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Abstract 

The dynamical origins of the EMC effect are studied. We conclude that a 
swelling in size of a bound nucleon as well as nuclear binding plays an impor¬ 
tant role in determining the parton distributions within a bound nucleon. We 
find that the distortion of nucleon structure functions in nuclei can be simply 
explained with a few fundamental nuclear parameters. 

PACS numbers: 13.60.Hb, 24.85.+p, 21.10.Dr 


The fact that the structure functions of bound and free nucleons are not equal is 
called the EMC effect [|lj . Recent accurate data Ml on nuclear structure functions 
impels us to reconsider the origin of the EMC effect. In this letter we report that 
this effect can be explained in a broad kinematical region using the idea of swelling 
and incorporating binding effects using only a few fundamental nuclear parameters 
but with a new understanding concerning these concepts. We expound on our ideas 
as follows. 

As the first step, we choose a set of parton distributions which can describe the 
structure functions of the free nucleon in the kinematic region of the EMC effect 
(where Q 1 2 ranges from less than 1 GeV 2 to a hundred GeV 2 ). One such model 
was proposed by Gliick, Reya, and Vogt ||. The model assumes that there exists 
a scale Q 2 = /a 2 which separates the perturbative regime ( Q 2 > /a 2 ) from the non- 
perturbative one ( Q 2 < A 2 ). All parton distributions are generated dynamically by 
evolution from a set of valence-like inputs (of the form Akr a P/v,g(^)( l—x)^) at /A. For 
example, in the leading order (LO) approximation with i 2 = 0.23 GeV 2 , the valence 
(u v , d v ), the total sea (5), and gluon (g) input distributions are given to be [|J, 


xu^(x) = 1.377x°- 549 P n , u (x)(1 — x) 3 ' 027 , 
xd^ (x) = 0.328x°' 366 Pn, d{x) (1 — x) 3 ' 744 , 
xS N (x) = 2x(u + d) = 2.40x°- 29 PAr ! 5(a;)(l — a;) 7 ' 88 , 
xgN(x,) = 35.8x 2 - 3 (l — a;) 4,0 , 
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where 


P N Jx) = 1 + 0.81v/i- 4.36a; + 19.4a; 3 / 2 , 

P N4 (x) = l + 1.14v/i +5.71x4- 16.9a: 3 / 2 , (16) 

Pn,s(x) = 1 + 0.31a; . 

The shape of the valence quark distributions in (1) is similar to that of the gluons, 
while the input sea quarks are much softer. This can be understood if we assume that 
the gluons are co-moving with the valence quarks whereas the sea quarks are mesonic 
at Q 2 = ii 2 , i.e., the input sea quarks correspond to the mesonic component of the 
nucleon at Q 2 < A 2 , so that their distribution will be softened by a convolution form 

!• 

We now discuss the nuclear binding effect at Q 2 = p 2 . Fermi motion is a global 
property of the bound nucleon. We regard the energy of Fermi motion as a part of 
the effective mass of the bound nucleon. This does not change the total momentum 
of the nucleon but only smears the parton distributions in the large x region. Since 
we are interested here primarily in the small and intermediate x region, we do not 
consider this effect in what follows. Therefore, the momentum of the struck nucleon 
in the lab frame is Pm = ( Mm — 26, p), where b is the binding energy per nucleon. 
The factor 2 is due to the fact that the nucleus with mass (A — 1) is required to be on 
mass-shell. Our approach is different from traditional models of the EMC effect 0 
since we do not need to compensate for the binding energy. Instead, we try by using 
the Weizsacker mass formula || to establish the connection between the binding effect 
and parton distributions in nuclei. According to this formula, the binding energy per 
nucleon arising strictly from the nuclear force is 

b = [1 - P s (A)}a vol + P S (A)^ = a vo \ - a sm A~ l/ * , (2) 

for A > 12 where P S {A) is the probability of finding a nucleon on the nuclear surface. 
We have ignored other contributions, especially the Coulombic one, which is not 
probed in DIS. Experimentally a vo i = 15.67 MeV, and a sur = 17.23 MeV. 

On the other hand, the attractive potential describing the nuclear force is from 
the exchange of scalar mesons. The energy required for this nuclear binding is taken 
away from individual nucleons which thus lose energy. We assume that the binding 
energy, b, corresponds to loss of energy of “mesonic” sea quarks in the nucleons. This 
means that the momentum fraction carried by the sea quarks in a nucleon bound in 
a nucleus at Q 2 = ii 2 will be reduced to 


S A (x,ii 2 ) 


K(A)S n (x, h 2 ) 

0 “ iiMSwV z )> 2 ) SN{x ’ t ‘' ) ' 


(3) 


Here (Sn )2 is the momentum fraction (second moment) of the sea quarks and we 
assume that the decrease in number of sea quarks due to the binding effect is propor¬ 
tional to the quark density. 

A bound nucleon will physically swell in size and this has been discussed in the so- 
called rescaling models, in which the rescaling of the input point, p 2 , is assumed [0. 
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However, we consider that the swelling of the nncleon only geometrically redistributes 
partons inside the nucleon and does not change either the value of this dynamical 
parameter or the existing number of partons at /i 2 . We try to describe the swelling 
effect using some universal principles as in The relative increase in the nucleon’s 


radius is 5a, where (Rn + A R(A))/Rn = 1 + 5a- Analogous to (2), we assume that 
the swelling of a nucleon on the nuclear surface is less than that of one in the interior; 
therefore, 

6a = [1~ Ps(A)}5 vol + P s (A)5 vol / 2 . (4) 


Here <5 vo i parametrises the swelling of the nucleon in the interior of a heavy nucleus and 
is a constant for nuclei with A > 12 and also for Helium (with P s = 1) since they have 
similar nuclear densities. Interestingly, the distortions of the density distributions and 
hence the changes in the three main parameters, N , a, and f3, in (1), due to swelling (if 
we assume P A} q(x) = P^ q (x)), can be simply determined by some universal principles. 
These cause the first three moments of the parton distributions in a free ((pv) and 
bound (q/f) nucleon (for q = valence, sea quarks, and gluons) at /i 2 to be related by 


(Qa(^ 2 ))i 

( Qa (^ 2 ))2 

(( Mh 2 )) 3 - ( Mh 2 ))^) 1/2 

(Mh 2 ))s - (u(h 2 ))l ) 1/2 


(Tv(T 2 )) 1 , 
(?v(/i 2 ))2 , 

1 + Sa ■ 


(5) 


The first two equations imply number and momentum conservation of partons and 
the last incorporates the swelling effect [|ItJ. Fig. 1 gives an example of the swelling 
effect for the parton densities in calcium. The momenta lost from the small and 
large x regions are transferred to the intermediate x region. The effect thus weakly 
enhances the distributions in the region 0.1 < x < 0.3 and results in “antishadow¬ 
ing”. Furthermore, this enhancement depends on the nuclear density and does not 
disappear at larger Q 2 . 

Finally, there is a further depletion of the sea densities at the time of scatter¬ 
ing. This is easiest to see in the Breit frame, where the exchanged virtual boson is 
completely space-like, so that the 3-momentum of the struck parton is flipped in the 
interaction. Hence, a struck parton carrying a fraction x of the nucleon’s momentum, 
Pjy, during the interaction time T mt = 1/u, will be localized longitudinally to within 
a potentially large distance A Z ~ 1/(2xPn)- Although the spatial extent of a single 
colored parton cannot exceed the range of QCD confinement, one can assume that a 
struck parton with x < Xq = 1/( 2 M 7 vdjv) can combine with a wee parton and form a 
colorless state T 0 with vacuum quantum numbers. Here d^r is the average correlation 
distance between two neighboring nucleons in the lab frame: 




P S {A) [R n + 2(i? ws 


(i - P.W) 


Rn 

~ir 


-Rat)] + 

+ 2(i?ws — Rn) , 


( 6 ) 


where Rn and i?ws are the nucleon and Wigner-Staiz radii. The observable value of 
the momentum of To is its average over the uncertainty time Tint and equals zero in 
the Breit frame. Therefore, Tq is a static scalar held. 
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Fig. 1 The effect of nucleon swelling on the calcium input distributions: the ratios of 
the modified to unmodified densities are shown for u v , d v , and S. 

We know that the correlation between two nucleons via scalar held exchange is 
the bound state problem. Therefore, in the DIS process, a nucleon can interact with 
other nucleons via T 0 in the long-distance. Because of this, a new binding effect will 
be experienced by the correlated nucleons, which we call the second binding effect. 

Consider a pair of nucleons which are correlated by To. The interaction causes 
the number of sea quarks S A (x, Q 2 ) in a nucleon to be further reduced to S A (x, Q 2 ). 
(Valence quarks are not depleted due to the requirement of quantum number conser¬ 
vation). Similar to (3), one can simply assume that the loss in number of sea quarks is 
proportional to the original density, i.e., S A (x,Q 2 ) — S A (x,Q 2 ) = /3S A (x,Q 2 ), where 
(3 is a constant. The total energy loss in a bound nucleon due to the second binding 
effect therefore is 

rx o 

U S (Q 2 ) = (3M n / xS A (x, Q 2 ) ~ (3M n (S a (Q 2 )) 2 . (7) 

Jo 

ffence, 

u,m u(p 2 ) ,, 

' M n (S a (Q 2 )) 2 M N (S K (^)h ’ 

U(p 2 ) = a vo i/6 being the binding energy between each pair of nucleons. 

Partons with momentum fraction x can overlap (n — 1) other nucleons, where 
n = l/(2M N d N x) = xq/x. The second binding effect arising from each overlap is 
simply additive due to the applicability of the superposition principle to the linear 
static held. The corresponding shadowing of the sea quarks is given by S A (x,Q 2 ) = 
K'(A)Sa(x,Q 2 ), where the depletion factor is, 

K'(A ) = 1, when x > xq] 

= 1 — 2f3(xox~ 1 — 1), when x A < x < xq, (9) 

= 1 — 2/3(xqx a 1 — 1), when x < x A , 
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where 2/3 = 0.037, xa = 1/(4:RaMn), and 2 Ra — IARa is the average thickness of 
the nucleus. We emphasize that the second binding effect acts on the intermediate 
state of the probe-target interaction and does not participate in the QCD evolution 
of the initial state. Since xq < 0.1, this is also a small x effect. 

In summary, the input nucleon structure function (for an isoscalar nucleus) is 
given by 


Fa{x,h 2 ) 



xu A (x, /U) + xd A (x, fl 2 ) 
+K(A)xS a (x, h 2 ) , 


( 10 ) 


where the average charge square factor is (e 2 ) = 5/18 and q A (x,/a 2 ) incorporates the 
effect of swelling on every input parton density, q N (x,/a 2 ). The structure function 
ratio of nucleons bound in two different nuclei, A and B , at Q 2 > /i 2 is given by, 


xu A (x, Q 2 ) + xd A (x, Q 2 ) + K'(A)xS a (x, Q 2 ) 
xu B (x, Q 2 ) + xd B (x, Q 2 ) + K'(B)xS b {x, Q 2 ) 


where Sa(x,Q 2 ) corresponds to the evolution of K(A)Sa(x, /i 2 
input sea density, depleted by the binding effect, to the scale Q 


i.e., 

2 


of the swelled 
Note that for 


Q 2 > /j 2 , the sea density also includes a finite contribution from strange quarks. The 
smearing effect of Fermi motion (at large x) is neglected in this work. 

In fig. 2 we give our predictions for the x, Q 2 , and A dependences of the ratios 
for He/D, C/D and Ca/D. The only free parameter to be fixed in computing these 
ratios is the value of <5 vo i. We find that <5 vo i = 0.15 best describes the various available 
data for these nuclei. 

We have also shown the ratio of the 6 Li to D structure functions in fig. 2. Unlike 
helium (which also has A < 12), lithium is a very loosely bound nucleus (whose 
density is less than half that of helium or carbon, while its radius is larger than that 
of carbon [§). Hence, 5a for Li can be expected to be much smaller than that for He 
($He = 0.15/2 according to (4)) and we take it to be 0.033. In all cases, the swelling 
effect dominates the intermediate x ratios while both the binding effects determine 
the small x ratios. 

The ratios R AB = F A /F B for A/B corresponding to C/Li, Ca/Li, and Ca/C are 
more sensitively dependent on the model of the EMC effect (as seen in fig. 3). We 
predict, in general, an enhancement of the ratio, R AB , at intermediate x which is a net 
effect of the swelling and the first binding effect, both of which depend on the nuclear 
density. On the other hand, the strong depletion at small x depends on both the 
nuclear radius and density. The weaker enhancement of the intermediate x ratios of 
Ca and C is because of their similar densities |TT[ as compared to the C/Li and Ca/Li 
ratios where the two nuclei have different densities. Our predictions are seen to be in 
excellent agreement with the NMC data [|5|, also shown for the sake of comparison. 
Note that we have computed these three ratios at the same values of (x, Q 2 ) as the 
data. We emphasize that, apart from the swelling parameter, <5 vo i, all parameters are 
fixed by a few fundamental nuclear parameters. 
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Fig. 2 The structure function ratios as functions of x for (a) He/D, (h) Li/D, (c) C/D, 
and (d) Ca/D. The dashed, full, broken, and long-dashed curves correspond to Q 2 = 
0.5, 1, 5, and 15 GeV 2 respectively. The data, i, m shown as open and solid circles, 
boxes and triangles correspond to Q 2 < 1, 1-5, 5-15, and > 15 GeV 2 respectively. 

The input (1) reliably predicts the free nucleon stucture function in the region 
Q 2 > 0.5 GeV 2 . Hence our dynamical model is also expected to predict correctly, 
not only the ratios, but also the values of the nuclear structure functions themselves 
in the same kinematic region. The detailed comparison between the theory and the 
data will published elsewhere. 

In conclusion we have shown that the EMC effect can be explained using swelling 
and binding effects but with a new understanding of these ideas. As the general 
properties of the nuclear force originated from the binding effect of nuclei in the 
history of nuclear physics, we expect that a new understanding of binding effects in 
the EMC effect will bring to light the nature of the nuclear force at the level of quarks 
and gluons. 
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Fig. 3 The structure function ratios for (a) C/Li, (h) Ca/Li, and (c) Ca/C. The 
solid curve is a smooth fit to our theoretical predictions at the same ( x, Q 2 ) as each 
available data point 1 The description of the data is the same as in fig. 2. 
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